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ABSTRACT 


A computer technique for obtaining near-field cross- 
section views of an ultrasonic field is presented. The 
technique is applied to ultrasonic transducer radiation 
patterns with sample contour drawings of several cross- 
sections. A spatial frequency spectrum formulation of 
diffraction theory is chosen for the computations which are 
performed with the aid of a fast Fourier transform routine. 

By experimentation an actual ultrasonic field was 
generated using a quartz transducer and the field was 
measured in a cross-sectional plane. A computer generated 
drawing of the measured data is compared to the drawing of 
a computer generated prediction of the field at the same 
plane. 

The computer model is applied to an array of transducers 
to investigate the resulting ultrasonic field, as a pre- 


liminary use of the technique for array design. 
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Ll. SEN ReDUCHION 


The work in this study was in connection with ultrasonic 
imaging systems studies at the Naval Postgraduate School. 
The study undertaken here was generated by the need for a 
large~area ultrasonic field as delineated in Ref. 1. One 
possible way of producing a large planar wavefront is by 
synthesis of a large array from smaller transducers. Due to 
the large distribution of variables involved in deSigning a 
near-field array, a facia practical means of producing 
near-field information was first required. A large body 
of literature is available for simple apertures in the near- 
tueidwand C. 0% MeKanney et alerecount the»history,of this 
research very concisely in Ref. 2. Exact expressions of the 
near-field are developed through essentially the same methods. 
They all involve finding a solution to the wave equation 
which satisfies certain assumed boundary conditions. The 
resulting integral equations are unwieldy even for the 
ee aperture shapes because the field cannot be 
expressed in terms of elementary functions. The ability to 
project planes of data in the near-field to other planes 
within the near-field is therefore a most valuable facility 
in terms of making maximum use of both theoretical and 
empirical data available for design purposes. Developing a 
model to make these plane-to-plane projections was the major 


task undertaken in this study. A starting point for planning 
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the work was a study of the imaging technique reported by 
meee BOyer et al [Ref. 3]. 

As a test of the model developed, the distribution of a 
Square aperture's diffraction pattern was predicted at 
various planes throughout the field. The model was then 
used to predict the field of a quartz transducer. The 
actual field distribution of the quartz transducer modeled 
on the computer was measured after setting up the necessary 
experimental apparatus. A piezoelectric probe attached to 
a mechanical scan system was used to collect the data. A 
preliminary use of the data was the comparison of a predicted 
and a measured field distribution in a plane within the near- 
field of the source. The data taken in this experiment 
should be useful for future work in this and related areas 
of imaging system studies. 

Time allowed the use of the computer me@e. to make pre- 
nies Fic of the near-field distribution of a transducer 
array. However, the use of the computer model for array 
design was done only in a preliminary manner by examining 
an ee array which was relatively simple to model. This 
allowed the method of design to be evaluated although a 
design problem itself could not be accomplished within the 


time alloted for this work. 





Ii. THEORY 


A. THE GENERAL SCHEME 

The purpose of this section is to present some of the 
fundamental concepts of an approach taken for computer 
modeling of the near-field of a transducer. It is not 
intended as a rigorous proof in the foundations of spatial 
frequency spectrum theory [Ref. 4], but rather a justifica- 
tion of the computer scheme used in this study. The overall 
scheme is one which allows the prediction of a complex near- 
field distribution in a plane on the basis of a known complex 
field distribution in any other plane perpendicular to the 
direction of propagation. Based on the assumption that the 
planar distribution of a sound field at the transducer takes 
on the geometry and motion of the transducer itself, predic- 
tions of the sound field distribution at various planes 
throughout the near-field were made with the computer model. 
Figure 2-1 is a schematic which is helpful in organizing the 
SPetatcicons involved in tne™sollowing subsections. ihe sub- 
sections contain an explanation of the symbols which appear 


in Page 2- 2. 


B. THE SPATIAL FREQUENCY SPECTRUM 
Let an isolated sound source be contained entirely within 


a source plane as shown below in Fig. 2-2. 








| 
Z=Z, plane Z=Z, plane 


Figure 2-1. The Propagation Model Sequence of Operations. 
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Figure 2-2. Two-dimensional Sound Source. 


Given that the source is a monochromatic pressure field in 
free space, and that x, y, and z are the cartesian coordi- 
nates of a typical field point, an expression of the pres- 


Sure distribution in the source plane is 


P(x,y,Z,t) = P(X,YrZ,) et? (Xr¥r Zo) emer (251) 


where t is time and w the angular frequency. Also defined 


are the following: 


= the complex pressure phasor 

P = the _—" pressure amplitude 

@ = the phase of the complex phasor 

Zo = the source plane position on the z axis 


For convenience an abbreviation in notion is used hereafter 


such that 


P_(x,y,t) = P@,y,z,t) (Zee ) 





where the subscript denotes the plane located at z. Assum- 
ing harmonic time dependence the expression for complex 


pressure is simplified to 

* Jo, (x,y) 

Po (x+y) = P_ (x,y) e (2:23) 
In viewing P (x,y) as a two dimensional distribution for a 
given z plane, the decomposition of the distribution into a 
spectrum of plane waves proves convenient. A decomposition 
of the distribution can be formed by use of a two-dimensional 


Fourier transform as given by Goodman [Ref. 5]. 
Br(fety) = f [B, Ory) expl-j2m (£xtfy) } dxdy (2.4) 


Goodman further shows that the original distribution P_ (x,y) 
can now be written in terms of its spatial frequency spectrum 


by use of the inverse Fourier transform as shown by 
Pee, y) - ff Pr(f,,£,) exp[j2m(f,xtfy)] af dt, (2.5) 


The idee che frequencies, f, and So can be related to a 
unique orientation of the propagation vector, k, for a given 
plane wave in 3-dimensional space. Where k, the propagation 
constant, is given by k = 2mT/A = 27nf/c, and f is the spatial 
frequency of a propagating plane wave in the direction of 
propagation; the vector, k, has magnitude, k, in the direc- 
tion of an outward normal from the planewave wavefront. Let 
the vector, k, be separated into its’ components in the coor- 
Gimease system by forming its dot product with a radius 


-> 
vector, r, such that 





k-r =k +k +k iil oa +f +f 236 
r= kx oe 22 = 2m(f x pl 72) (2-6)) 


Eq. 2.6 can be written in terms of the x-y plane spatial 


frequencies alone by making the following substitution: 


Mm the form given iin Eq. 2.5, P (x,y) expresses the distri- 
bution of sound pressure over the observed plane, z, as a 


superposition of plane waves and evanescent waves. Plane 


waves are implied when a = a “s + and evanescent waves are 
r 
. 2 2 I 
implied when f° + f° > —> [Ref. 6]. 
x y 12 


C. PROPAGATION OF THE SPATIAL FREQUENCY SPECTRUM 

The propagating waveform must satisfy the homogeneous 
wave equation in 3-dimensional space. P_ (x,y) in the frequen- 
cy spectrum formulation as given in Eq. 2.5 is required to 


satisfy the Helmholtz equation 
VP. +k P_ = 0 2.7) 


The solution of the Helmholtz equation leads to a second 
order differential equation in the variable z. Its solution 


1s given in Ref. 5 as 


a P co ip 1 2 | 
P (X,Y) Lf Pp ey ty exe pS fe - exe ances ace, 


(2.8) 
Comparing Eq. 2.5 and Eq. 2.8 it is apparent that Eq. 2.8 


can be considered the two-dimensional Fourier transform of 


that portion of the expression under the integral gven by 


10 





_ : il 2 2 
' — ———_ oF = 
Potty ty) exp sare 2 t.. t, | 
This expression suggests that propagation could be modeled 
in the Fourier transform domain using as a transfer function 


: . ee eee 
ee = exp sams eo fi . | (2.9) 


by merely taking the inverse transform of Pe xt y) Bt rt) 
where z in Eq. 2.9 is taken as the distance between the 


reference plane and the observed plane. 


D. TWO-DIMENSIONAL LINEAR SYSTEM 

It remains to be shown that the propagation from one 
plane to another can be considered an overall linear system 
to justify the transfer function approach to modeling. In 
a linear system the total response is equal to the combina- 
tion of elementary responses, or mathematically stated the 
Superposition property holds [Ref. 7]. In order to demon- 
strate this property of frequerryspectrum propagation, it is 
convenient to express the waveform at an observed plane 
Z = Z, in terms of the waveform at a reference plane, Z=Z.. 

To formulate a.wave distribution at the z = Zo plane as 
shown in'Fig. 2-3, a substitution of f., into Eq. 2.8 and a 
regrouping of that expression leads to a planar distribution 
given by 

00 

Be. (eeny 2) ma Py (Ege fyexp [2m (fx tty otf z)ldf,af, (2.10) 
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Figure 2-3. Boundaries and Distributions of a 
Propagated Wave Distribution at 
Two Planes. 


When the Fourier transform is denoted by F{ } and the 


inverse Fourier transform by eee }, Eq. 2.10 becomes 


P (27.40) 


3 : o j2nt 2. 
P (x_,y_.) = F -(P' (££) e i 
a ero Zn Ox Y 


Taking the Fourier transform of Eq. 2.11 and reordering the 
equation gives 


a eo i exp[j2m£, 24] = r{P, (or¥o)| (272) 


which can be equated to solve for P) (iets) as given by 
O 


| = —4 D 
Be (exe fy) = exp[-j2rfz.] [J Pz (Xor¥o) 


x exp ([-j 2m(f,x +f y,)] dxdy (2543) 


At the plane z = z, write the expression for P_ (x3yys) by 
d Z 4 are 


12 





Swestituting Eq. 2.13 for = asta im eg. 2ese The result- 


ime expression foreP x.,y.) in termseof P_ (x,y) is: 
Z 4 a= d Zo 


Baya “ff exp[-j2mf, (24-2,) JJ Pe oY! 


x exp [-j2m (fx +f y,)] 


x-exp[j2m(fx5+f ee oe ee (2.14) 


y 


By interchanging the order of integration of this equation 
and by making the simplification in notation that 


ise] 


K 59 (XQ XY Q7 YQ) = {f exp | j2m[£, (x4-X4) 


+ £ (yg) + £,(zg-2,) }} dt ,ae, C2". L5) 
the final form of the expression of Ee at the z = Z 3 plane 
d 
in terms of at the z = Z 6 plane can be expressed as 


P24 “ava! = JP Pe Xoo) Rao Xa" XorVa"¥o) MX gA¥g (2216) 


Systems whose output can be expressed in this form as a 
function of the input are recognized to be linear systems 


[Ref. 8]. Thus the wave distribution at the z = plane 


8 
is given in terms of the distribution at the z = Zo plane by 
Eq. 2.16 when the system transfer function is given by 

Bey. 22>. - 


In the computer model developed for Eq. 2.9 evanescent 


waves were neglected and the transfer function was therefore 


13 





bandlimited to those frequencies within a circle of radius 
1/rX. This at first appears to limit the class of pressure 
distributions which can be considered to those which have 
bandlimited spatial frequency distributions. However, evan- 
escent waves decay exponentially with increased distance in 
the z-direction, and their effect on the field is insignifi- 
cant .i1n-planes braaelle ot ee a few wavelengths from the source. 
The region of the field where the evanescent waves have a 
Significant effect on the pressure distribution is sometimes 
called the very-near-field. Acoustic imaging systems 
usually reconstruct the image from the pressure distribution. 
Reconstruction of an image requires that the evanescent 
waves be accounted for in the superposition integral of the 
type given in Eq. 2.16. Since the evanescent waves do not 
propagate, the near-field and far-field distributions are 
independent of their presence. Neglecting the evanescent 
waves, therefore, did not restrict the ability to make 
plane-to-plane predictions in the near-field as long as the 
reference distribution was measured at a distance greater 
than a few wavelengths from the transducer. Expressing 

Eq. 2.9 in a more appropriate form for the model which 


neglects the evanescent waves 


; | 1 2 ), 2 C xX. 
exp j271z [52 t.. = swat fy. + ty < 52 | 
H(f,-f.) = (2.17) 


0 . : otherwise 
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Be WERE DISCRETE FOURIER TRANSFORM 

Let X and Y be the sample intervals in the spatial x and 
y coordinate directions respectively. Let QS and Ry be the 
corresponding sample intervals in the spatial frequency 
plane. The discrete Fourier transform (DFT) over a region 
N and M samples long in the x and y directions respectively 


is defined as 


N-1l M-1l 
F(KQL La a pet (oe ey) 
Y n=0 m=0 
x exp 7 (neha tme ea) (Zero) 
where Deis < N=1 and D Saae SaeMi= 1 
and 
_ 2 _ 2m 
an » Cy Tee 


These equations are in the form of those which were imple- 
mented with the fast Fourier transform: (FFT) algorithm 
called "FOUR2" [Ref. 9]. The Fortran statements of these 
programs are included in Appendix A. The FOUR2 subroutine 
was used to perform the transformation to and from the 


Spatial frequency spectrum domain. 


F. PROPAGATION BY TRANSFER FUNCTION 


When Eq. 2.17 is placed into a discrete format as 


follows: 
: 1 2 2 ; 
exp |j27z I = (KQ,) - DEF) : 
eo. Qe) = (en)? + (29.92 <« + + (2.19) 
x’ VY x VY x2 


0 : elsewhere 
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or more simply written 


; il Zz oe 2 ; Boe Z 
exp [j21z 7) - — (kK°+2°) af (k™+2 De <I 
H = (2.20) 
im % 2,2 2 
0 safe +2 2X @-- 1 
XY a 


The propagation algorithm as given in Appendix A was used to 
perform the multiplication of the transformed input data by 


the transfer function (Eq. 2.20). 


G. ANALOG-TO-DIGITAL CONVERSION 

When predicted patterns of the waveform at various planes 
were calculated on the basis of their aperture function, the 
data was arranged in the shape of the aperture with uniform 
amplitude and phase samples everywhere within the aperture. 
However, the anGMe~aatva collected in the experimental por- 
tion of this analysis had to be converted to sampled data 
before being utilized in the propagation program. The 
analog-to-digital, A/D, conversion process was done by the 
Simultaneous utilization of an XDS 9300 digital computer and 
a COMCOR Ci 5000 analog computer. A lengthy description of 
the facilities and equipment used, the A/D conversion pro- 
cess, and the seven-to-nine track conversion was a necessary 
step in making the digital data of the XDS 9300 computer 


compatable with the IBM 360/67 computer requirements. 


H. AN INTERPOLATION OF DIGITAL DATA DISPLAY 
The subroutine CONTUR [Ref. 11] included as part of 


Appendix B, was used to produce separate analog displays of 
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the amplitude and phase at a desired plane. In this routine 
a pseudo three dimensional effect is presented by contour 
lines. Coarse analysis of an entire cross-section of the 


field can be performed from contour drawings. 


I. FLOWCHART OF STEPS IN COMPUTER OPERATIONS 
A flowchart of the computer operations in the order 
performed and a description of the equipment used in con- 


junction with each step is given in Fig. 2-4. 


oi. INHERENT ADVANTAGES AND DISADVANTAGES OF THE DFT AND FFT 

Since Cooley and Tukey first reported the algorithm 
which they developed for the FFT [Ref. 12], the amount of 
literature which has been published on applications of this 
technigue of computation and its advantages and disadvantages 
would fill volumes. Much insight into the FFT technique can 
be gained from Refs. 13-15. Much less has been written 
about applying this technique for processing of data in the 
two dimensional spatial frequency domain. Reference 13 was 
the most germane discussion found for this study. It 
discusses the problems commonly called leakage and aliasing 
in articles which treat time-frequency analysis problems 
(Ref. 14). 

The problem of leakage can Senetanas work as an advan- 
tage in spatial frequency analysis. Leakage is an approxi- 
mation error which arises from a finite record length 
representation of an infinite length raged. Its affect on 
the desired spectrum is equivalent to an undesired convolu- 


tion of a data window with the true frequency spectrum. The 


IL 
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Figure 2-4. Flowchart of Computer Operations for 
Propagation Modeling. 
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window occurs because of the geometry of the termination of 
the record. In this study a transducer array was in reality 
a truncated source and propagation of the data window 
resulted in a diffraction effect of the window geometry. 
Diffraction was modeled using this effect when the square 
shape of the sampled lattice was suitable. When the square 
window did not adequately model the geometry desired, the 
simplest method of modeling the edge diffraction was found 

to be one of leaving a border of zeroed data around the data 
field to be transformed. This method reduced the diffraction 
effects of the square sample lattice at the cost of a greatly 
increased memory requirement to model a given transducer. A 
discussion of the truncation approach to handling various 
geometries iS given in Ref. 17. 

Aliasing is a problem associated with modeling a con- 
tinuous Fourier transform (CFT) by a discrete Fourier trans- 
form (DFT). It is, therefore, an inherent disadvantage of 
the DFT. The term “aliasing" refers to high frequency 
Semponents of the spatial function which are incorrectly 
represented due to widely spaced sampling. The effect is 
worse than a loss of the high frequency components since 
it results in erroneous coefficients in the lower fre- 
quencies of the spectrum. The effect can be avoided by 
sampling the waveform with the optimum spacing [Ref. 18]. 
The sample spacing when determined by the two-dimensional 
sample theorem frequently led to excessive numbers of 


samples necessary to represent a given transducer. Since 
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computations were computer memory limited, a corresponding 
limit was placed on the total number of samples available 
to represent a cross-section of the field. A compromise in 
choosing sample spacing and border width around the repre- 
sentative data was necessary. 

The disadvantages of the FFT associated with modeling 
the diffraction equations were insignificant when compared 
to those of the DFT. That is to say that the errors in 
modeling of the diffraction came from making the transition 
Erem aeCFT to a DFT not im the FFT algorithm itself. The 
worth of the algorithm is measured in terms of speed, 
efficiency, memory required and round-off error. Of these, 
the requirement for a large memory is the only significant 
disadvantage. Although this study was seriously limited by 
the computer memory available, there are means to circumvent 
this problem by using more computation time. This is a 
matter of optimizing the algorithm to the particular Baca 
of the propagating waveform that is of interest. For this 
work, three versions of the algorithm listed in Appendix B 
were required. The subroutine "FOUR2" makes very effective 
use of time in performing the calculations of the FFT, how- 
ever, the total number of sample points (record length) must 
be a power of 2. Reference 19 gives an algorithm which 
can transform many different record lengths at a moderate 
cost in increased computation time. Implementation of this 
"mixed radix" form of algorithm is strongly recommended for 


future studies of transducer sound field analysis. 
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Til. BXPERIMENTAL INVESTIGATION 


Experimental data was taken for two primary reasons. 
First was the desire for a means of evaluating the computer 
predictions. Second was the desire for actual seta vice the 
assumed data at the face of the transducer as mentioned in 
Section II to evaluate the assumptions of the model itself. 
A mechanical scan probing system was constructed to measure 


the complex sound field of a transducer to be described. 


A. THE WATER TANK 

The water tank used in this experiment was constructed 
of plywood and lined with a one inch thickness of poly- 
styrene. The inner dimensions provided for a column of 
water 4% inches deep, 4 inches wide and 45 inches long. A 
cross section of this column was not large enough to con- 
sider the propagating wave entirely unbounded. However, it 
was found by experiment that the difference in nearby rigid 
and nearby release surface measurements were not significant 
until the hydrophone probe was within approximately ten 
wavelengths of the boundary. All data was taken with this 
consideration limiting the field of data. The length of the 
tank used was considered adequate to approximate an infinite 
medium at the frequency used. The tank was lined with poly- 
styrene which was a composition of pliable thermoplastic 
with a high absorption coefficient. Its impedance is very 


@lose to that of air for Which p.c, = 415 rayls. 


2a 


B. “HE HYDROPHONE PROBE 

In design of the hydrophone (Fig. 3-7a), an attempt was 
made to minimize reflections affecting the free field 
measurement. The active face of the probe is 0.04 inches in 
diameter (0.83 wavelengths @ 1012 kHz). A thickness vibrator 
type receiver of PZT-5 ceramic composition is produced by 
Valpy-Fisher. .It is encapsuled with a bronze impedance 
matched backing and provided with a shielded electrical 
connector by the manufacturer. This type probe was mounted 
for use aS a hydrophone in a plexiglass holder and coated 


with neopreme. 


Cc. THE TRANSDUCER 

The source transducer is shown in Fig. 3-lb. The crystal 
Shown is a 2 inch by 2 inch square transducer of X-cut quartz. 
This transducer is designed as a thickness vibrator polished 
fr a fundamental frequency of 1 MHz with good overtone opera- 
tion characteristics. The peak output signal from the probe 
was obtained when driving the transducer at 1012 kHz. Con- 
sequently that frequency was chosen for all data collected. 
The small offset from the transducer's fundamental resonance 
was necessary due to the mis-~match in resonant frequencies 


of transducer and probe. 


D. THE MECHANICAL SCAN EQUIPMENT 
The mechanical scan equipment mounted above the measure- 
ment tank is shown in Fig, 3-2. The scanning system became 


an important concern in fabricating the measurement apparatus 
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because three rigid system requirements depended upon the 
scan mechanism. The first was the requirement to insure 
that phase information be recorded accurately to within +10 
degrees. This meant that the recording plane would have to 
be flat to within 0.04 millimeters at the frequency used of 
1012 kHz. A similar precaution was taken by Boyer et al as 
noted in Ref. 3. The second requirement was that the scan 
system be equipped to measure position accurately to within 
a small fraction of a wavelength. Still a third requirement 
was that it provide a means of taking data rapidly enough to 
collect a large quantity of data in a reasonable length of 
time. The heart of the scan system was a mechanical scan 
unit taken from a disassembled fire control computer. This 
unit, driven by a variable speed motor, positioned the probe 
horizontally across the width of the tank. A manually 
operated slide assembly positioned the probe in the vertical 
dimension to vary its depth. The coordinate position of the 
probe within the cross-sectional plane was calculated by the 
advance per turn of the threads on the positioning mechanisms. 
The ress of turns was recorded by hand from a mechanical 
counter in the vertical, and by a potentiometer in the hori- 
zontal direction. Analog records were taken for each depth 
as the motor driven assembly scanned the width of the tank. 
Uniform depth increments were used between records so that 
ultimately a square matrix of equally spaced data points 
would result. The assembled system produced data which shows 


a phase variation across the face of the transducer of 20 


25 





degrees peak to peak. The phase information recorded in 
Fig. 3-3 was taken with the probe moving horizontally across 
the width of the transducer 9.1 wavelengths in front of it, 
at a point midway between the upper and lower edges of the 
transducer. There is an overall shift in the mean phase 
from the left edge to right edge of the transducer of 
approximately 30.degrees. This indicates a misalignment of 
the transducer and recording planes by 0.12 millimeters 
maximum. Phase records were used in this manner to align 
the transducer to the recording plane in both the horizontal 
and vertical aspect. The maximum deviation in phase about 
the mean was +10 degrees when taken directly in front of the 
transducer. On the basis of this and similar records it was 
estimated that phase measurements were accurate to within 
+10 degrees at the frequency used of 1012 kHz. Turns of the 
scan positioning mechanisms were recorded accurately to 


within 1/36 of a turn with a resultant position accuracy On 
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Figure 3-3. Phase Variation across the Face of a Quartz 
Transducer Radiating at 1012 kHz. 
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approximately 1/50 of a wavelength. The sample spacing was 


9/10 of a wavelength at uniformly spaced intervals. 


E. ELECTRONIC MEASURING EQUIPMENT 

The general layout of electronic measurement equipment 
is shown in Fig. 3-4 and a photograph of the physical layout 
is given in Fig. 3-5. The electronic equipment layout was 
an assembly of standard pieces of equipment. A 
general description and specification summary is given 


in Appendix B for the major components. 


F. EXPERIMENTAL PROCEDURE 

It was pertinent to the data reduction process that the 
diameter of the hydrophone face was less than the sample 
spacing, thus avoiding overlapping of the samples. Had 
imaging been the desired use of the data collected, it would 
have been ideal to sample the field at half-wavelength 
intervals. The scanning system described in this experiment 
could have been used for imaging if a more complex data 
reduction process were used [Ref. 3]. The distance between 
the transducer and the recording plane was adjusted by 
moving the transducer which was mounted on a traveling 
holder. Separation of transducer and probe was measured by 
a scale graduated in 1/64 inch. If time had allowed repeated 
experiments of this type, a mechanically driven threaded 
shaft and traveler system would have been used for convenience 
in observing ultrasonic wave behavior in the direction of 


a 


propagation, 
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‘Figure 3-4. Electronic Equipment Layout. 
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G. TYPICAL EXPERIMENTAL RESULTS 

Pagure 3-6 sh@ws six typical records as recorded at 
successive depths by the analog type scanning equipment. 
Amplitude, phase, and potentiometer outputs are recorded 
Simultaneously on separate channels. The experiment was 
designed to generate data for a 64 by 64 sample lattice, so 
the potentiometer output has a spike corresponding to each 
of 64 equally spaced positions of the hydrophone along the 
width of the record. This channel was used in the A/D 
conversion to trigger the sample circuit. Figure 3-7 
shows the analog playback of the sampled data as recorded 
on magnetic tape after the A/D conversion. Comparison of 
the measured record and the converted record shows good 


agreement. 
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Six Typical Records of the Sampled Data Records Showing Amplitude 
and Phase (reading top to bottom) as Played Back from a Magnetic 


Tape in Analog Presentation. 


Figure 3-7. 





IV. EXPERIMENTAL vs THEORETICAL RESULTS 


A test of the computer program was the generation of a 
far-field diffraction pattern. This was done by using a 
small aperture whose diffraction pattern could be considered 
to be in the far-field at distances well within the near- 
field of the sample lattice used. The effect of the sample 
lattice was introduced in the preceding section in the 
discussion of the DFT where the scheme of making the perture 
to be modeled much smaller than the "data window" was pro- 
posed. This scheme allows increased latitude in the range 
of z for which the diffraction model can be considered valid. 
If a sufficient number of discrete samples are taken to 
accurately represent a square aperture, and a large border 
of zeroed data is left around the aperture, then ne window 
does not disturb the desired diffraction pattern. It was 
convenient to consider this effect in the following light: 

The very-near-field of an aperture is essentially a 
parallel projection of the aperture. The diffraction of a 
small aperture diverges very rapidly and the far-field is 
reached rather quickly in moving in the z-direction. The 
diffraction pattern of a large aperture diverges more slowly 
and the far-field is reached at much greater distances of 
propagation. (The same radiation frequency for both aper- 
tures is assumed.) The terminology "near-field of the data 


window" will be used to refer to the region in which the 
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modeled aperture is valid. (Equation 2.9 which is the exact 
and continuous form of the propagation model has no such 
limitation on the region of the z axis in which it is valid. 
The limitation is a result of a finite sample lattice and 
with a large enough sample space the diffraction of any 
aperture can be projected into the far-field with this 
model. ) 

W. K. Pratt and H. C. Andrews [Ref. 20] have generated 
computer far-field diffraction patterns which were used for 
comparison to test the computer program for this work. 
Figure 4-1 shows a typical diffraction pattern of a single 
Square aperture. The amplitude variations of the field 
(ranging from 0.0 to 3.1 units) are present in five contour 
levels, for a plane 96 wavelengths (964A) from the aperture. 
The divisions between levels are 0.7 units apart. The 
width of the aperture, a, wee 6X and the width of the data 
field was 38 across. The pattern shown was as expected 
from the analytic solution for a square aperture [Ref. 5]. 

A variety of criteria are found in the literature for deter- 
Ming the transition from the near-field to the far-field 
region. According to Goodman [Ref. 5], this transition can 


be given as 


jh 
ke x FY) 
Zz 5 (4.1) 


where z and k are as previously defined and x and y are 
measured from the center of the aperture as shown in figs. 
2.2 and 2.3. Writing Eq. 4.1 in terms of sampled data for 


a square lattice gives the criterion 
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Z>> ==> NXA (4.2) 


From this equation the maximum z for which the window does 
not appreciably interfere with the observed diffraction 


pattern is taken as 
Z = =~ N'A (4.3) 


For example, the pattern shown in Fig. 4.1 is a square 
lattice of 32 samples giving a far-field for the data window 
beginning at z = 1340. Applying the same equation (Eq. 4.1) 


to the maximum x and y, where x = 


a ae = N/16, the far- 


field of the aperture is where 


2 


ie a) or Zz >>17.4A 


21 
me ‘Te 


mmvcds expected, therefore, that Fig. 4-l would be a reason- 
ably accurate ee ee OP REiordicat lewd Palawan Ox ict 
square aperture out to z = 1340). 

In the preceding discussion the word aperture was used 
instead of the word transducer because of the type of distri- 
bution assumed at the transducer plane. When the following 
assumptions are made the transducer diffraction pattern is 
based entirely on its aperture function: 

1. The motion of the transducer has uniform amplitude 
over its entire face. 

2. The motion of the transducer is in phase over the 
entire face. 

3. The motion of the medium immediately adjacent to the 


face of the transducer has the same motion as the transducer 


itself. 
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4. The medium outside of the edges of the transducer 


does not move due to the motion of the transducer. 


The program was tested by its generation of aperture diffrac- 
tion patterns. The ability to predict transducer diffraction 
patterns on the basis of an aperture diffraction pattern will 
be discussed in light of the experimental data. 

The transducer shown in Fig. 3-lb-was modeled on the 
IBM 360 computer and the field at z = 73A was predicted 
based upon its square aperture function. Figures 4-2 and 
4-3 are the predictions of amplitude and phase eect e 
Note that the transducer was a 2 inch square sheet of quartz, 
while the active face assumed in the model was a 1 3/4 inches 
square aperture. This Wasea Gerrection for that portion of 
the transducer obscured by the mount. Additionally, the 
restriction in motion offered by the mount at the edges 
violates the assumption of uniform amplitude. This effect 
can be seen from the amplitude record made at 9.1A in front 
of the transducer. Figure 4-4 shows a typical amplitude 
measurement at z = 73A. The assumed shape at the trans- 
ducer is shown by the heavy dotted line superimposed over 
enevdata record. Figures 4-5 and 4-6 are contour. plots of 
the sampled data at z = 9.1\ as recorded after the A/D con- 
version. The disparities in the predicted and measured 
contour mtterns is attributable to the following inaccuracies: 

1. The wavefront and measurement planes were slightly 
misaligned as seen from the phase record of Fig. 3-3. In 


order for the fourfold symmetry of Figs. 4-2 and 4-3 to be 
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Figure 4-2. Predicted Near-field Diffraction Pattern 
(amplitude) of a Square Transducer. 
a7rxe= 30, 2 = TSA. 
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Figure 4-3. Predicted Near-field Diffraction Pattern 
(phase) of a Square Transducer. 
a/AX = 30, z2 = 73d. 
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Figure 4-5. Measured Near-field Diffraction Pattern 
(amplitude) of a Square Transducer a/dA = 30, 
Z= 73d. 


41 








Figure 4-6. Measured Near-field Diffraction Pattern (phase) 
of a Square Transducer a/dA = 30, z = 732. 
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present in the measured pattern the source must be aligned 
symmetrically with respect to the measurement plane. The 
misalignment also caused the Beretta waveform to shift 
its center with respect to the center of the sample lattice. 

2. Data were not taken of the entire plane. To con- 
serve time in the recording and processing of data, only half 
of the plane was sampled; symmetry was assumed and the half- 
plane data was folded about a horizontal through the center. 
This gave a false symmetry to the data in view of the aiige= 
ment difficulty. 

3. Recording errors in the measurements taken were on 
the order of +10 degrees for phase and +4% for amplitude, 
due to inaccuracies of the recording equipment. The A/D 
conversion process introduced additional error, but no 
information was available on the accuracy of the program 
used other than the visual comparison from displays such 
as the one shown in Fig. 3-7. 

4. Inaccuracies in the free-field assumption caused 
by the probe hydrophone, and by the boundaries in the tank 
were neglected. 

5. Inaccuracies were inherent in evaluating a trans- 
ducer's diffraction pattern based on its aperture function 
alone. 

It is impossible to evaluate the validity of modeling 
the transducer by its aperture function on the basis of the 
comparisons presented here. If the data reaerawa at 9.1A 


were propagated to the plane at 73A and compared to the 





measurements made there, the composite effect of the inaccu- 
racies listed in (1) to (4) above could be determined. It 
was concluded for the purpose of this experiment that the 
aperture function model was sufficient for moceieng emat] 
transducers, and was found to be an adequate model when the 
ratio of a/A was less than 10. 

Many acoustic imaging systems of interest require planar 
ultrasonic fields at frequencies of 1 MHz and higher. An 
array design approach to supplying this need was the motiva- 
enon for developing the computer prediction capability. As 
a preliminary study to the design problem a 3 x 3 array of 
Square transducers was chosen. After some experimenting with 
the program it was found that at least 25 samples per trans- 
ducer was required to produce an acceptable far-field 
diffraction pattern for a square aperture. Each square 
transducer was represented by an aperture 6 wavelengths on 
a side. The spacing between these transducers had to be 
wide enough so that interference between adjacent transducers 
did not occur until a distance where the far-field pattern 
of the individual transducers had already developed. This 
effect is evident in the pattern shown in Fig. 4-7, where 
the sound field can be described as a sum of nine spatially 
separated far-field patterns. When the array was propagated 
to greater and greater distances, the far-field patterns | 
expanded and merged together in a way which made the indi- 
Weelal faretield patterngdistinguischable ian the overall 


field pattern. When the array was propagated to that distance 
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Figure 4-7. Diffraction Pattern for a transducer Array of 
9 Uniformly Spaced Elements. The elements are 
Squares,6 wavelengths on a side. The pattern 
is a prediction for a plane 60 wavelengths 


from the source plane. The contour levels are 
2.8 Units apart. 
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Figure 4-8. A Predicted Diffraction Pattern for a Transducer 
Array of 9 Uniformly Spaced Elements. he 
figure is one quadrant of a pattern which has 

four-fold symmetry. The elements are squares 
6 wavelengths on a side. The pattern resents a 
plane 600 wavelengths from the source plane. — 
The contour levels are 0.28 units apart. 
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where adjacent elements began to interfere, the field 
diffraction pattern became very complex and all similarity 
to the far-field pattern of a square was lost. Instead 

the pattern looked somewhat like the near-field diffraction 
pattern of a finely sampled square aperture. Figure 4-8 

is an example of this effect which the author has called 
"“the~near-field diffraction of far-field sources". MThis 
nomenclature alludes to the fact that the individual trans- 
ducers were in the far-field according to the criteria of 
Eq. 4.1 while the outer dimensions of the composite array 
placed the source as a whole in the near-field. This figure 
is the result of contouring only one quadrant and using a 

128 x 128 sample lattice. (Note the symmetry of Fig. 4-8, 
and the improved resolution gained by the larger number of 
samples used.) 

The use of the computer model appears to have much 
promise for designing an array to produce a desired waveform. 
There are a large number of variables which determine the 
form of the array diffraction. It was found, however, that 
the —" of trials which must be run can be significantly 
narrowed by defining dimensionless qiantities which model 
many different combinations of element size, element spacing 
and element shape. Variables which were not looked into in 
this preliminary study were the number of elements in the 
array, and the possible variations in driving the individual 
elements in various phase and amplitude combinations. If 


all feasible combinations of the variables mentioned were 
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to be examined for a design, it would be exhaustively time 
consuming and expensive. A major advantage of this approach 
is the intuitive feeling the user gains by running a few 
simple cases using well known aperture shapes. Familiarity 
with the geometry of some two-dimensional Fourier trans- 
forms [Ref. 21] coupled with a feeling for the general trends 
of propagation in the spatial frequency domain, would allow 
the designer to eliminate a large class of possible variables 
and to concentrate on those which have some intuitive promise. 
Facilities are available at the Naval Postgraduate School 

for immediate graphical response to a change in variables 
using the graphics display terminal which presents data on 

a rectangular cathode-ray tube. This type of data presenta- 
tion used with a routine similar to the CONTUR subroutine 
would make the large number of variables involved manageable. 
The IBM/360 version of the propagation program as presented 
in Appendix B requires 250K core storage with the CONTUR 
subroutine included. ThisS version requires approximately 

56 seconds of model 67 cpu time to compile and calculate 

the data for the first projection of a square matrix of 

64 x 64 samples. Each additional propagation of the data 

to another plane takes approximately 20 seconds. The time- 
memory requirement of the program appears to be adaptable to 


use with the graphic display units. 
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V. CONCLUSIONS 


A recently developed technique for imaging by reproduc- 
tion from object diffraction patterns has been applied to 
analyzing the sound source itself. The technique is not 
limited to ultrasound, but high frequency imaging applica- 
tions were of particular interest for this study. Conse- 
quently, the experimental measurements were conducted at 
1 MHz. Although most of the concepts and arguments can be 
applied to frequencies well below 1MHz, such applications 
were not considered herein. 

A technique of computer modeling of the diffraction 
equations has been emphasized. The capability of the 
computer facilities at the Naval Postgraduate School to 
perform such diffraction calculations in a fast concise 
format has been accented. 

The use of this technique in designing an array was a 
driving incentive in the work reported, and from the results 
of the preliminary array model examined it appears to be a 
worthwhile approach to production of a planar field of 
ultrasound. 

The experimental portion of this work produced data for 
an elementary comparison with a computer prediction. The 
data has many applications for future work to give infor- 
mation about the transducer used, about the capability of 


the computer program, and about imaging system capabilities. 
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APPENDIX A 
PROPAGATION PROGRAM — IBM 360/67 COMPUTER 
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GENERAL DESCRIPTION AND SPECIFICATION SUMMARY 
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